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Based upon the rate equations for the photon distribution function obtained in the previous paper, 
we study the formal solutions in three different representation forms for the Sunyaev-Zeldovich 
effect. By expanding the formal solution in the operator representation in powers of both the 
derivative operator and electron velocity, we derive a formal solution that is equivalent to the Fokker- 
Planck expansion approximation. We extend the present formalism to the kinematical Sunyaev- 
Zeldovich effect. The properties of the frequency redistribution functions are studied. We find that 
the kinematical Sunyaev-Zeldovich effect is described by the redistribution function related to the 
electron pressure. We also solve the rate equations numerically. We obtain the exact numerical 
solutions, which include the full-order terms in powers of the optical depth. 
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I. INTRODUCTION 

The Sunyaev-Zeldovich (SZ) effect [U, i, i, Q, which 
arises from the Compton scattering of the cosmic mi- 
crowave background (CMB) photons by hot electrons 
in clusters of galaxies (CG), provides a useful method 
for studies of cosmology. For the reviews, for example, 
see BirkinshawP and Carlstrom, Holder, and ReeseQ. 
The original SZ formula has been derived from the Kom- 
paneets equation|7| in the nonrelativistic approxima- 
tion. However, recent x-ray observations (for example, 
Schmidt et al.Q and Allen et al.0) have revealed the 
existence of high-temperature CG such as fcsTe ~20keV. 
Wright [13] and Rephaeli and his collaborator [ll|,[l2| have 
done pioneering work including the relativistic correc- 
tions to the SZ effect for the CG. 

In the last ten years remarkable progress has been 
made in theoretical studies of the relativistic corrections 
to the SZ effects for the CG. Stebbins[l^ generalized the 
Kompaneets equation. Challinor and Lasenbvp3| and 
Itoh, Kohyama, and Nozawafisj have adopted a rela- 
tivistically covariant formalism to describe the Compton 
scattering process and have obtained higher-order rela- 
tivistic corrections to the thermal SZ effect in the form 
of the Fokker-Planck approximation. Nozawa, Itoh, and 
Kohyama^l^] have extended their method to the case 
where the CG is moving with a peculiar velocity with 
respect to the CMB and have obtained the relativistic 
corrections to the kinematical SZ effect. Their results 
were confirmed by Challinor and Lasenbv[l7| and also 
by Sazonov and Sunvaev|18i [Toj . Itoh, Nozawa, and 
Kohyama [20j have also applied the covariant formalism 
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to the polarization SZ effect^, |4]. Itoh and his collabo- 
rators (including the present authors) have done exten- 
sive studies on the SZ effects, which include the dou- 
ble scattering effect [2ll|. the effect of the motion of the 
observer [22]. high precision analytic fitting formulae to 
the direct numerical integrations [H, and high pre- 
cision calculations [IE, [IB]- The importance of the rela- 
tivistic corrections is also exemplified through the possi- 
bility of directly measuring the cluster temperature using 
purely the SZ effect [27*1. 

On the other hand, the SZ effect in the CG has been 
studied also for the nonthermal distributions by several 
groups [H, [2^, [s^l- The nonthermal distribution func- 
tions, for example, the power-law distributions, have a 
long tail in high electron energy regions. Therefore the 
relativistic corrections for the SZ effect could be more 
important than the thermal distribution. 

The relativistic SZ effect has been studied in several 
different approaches. Wright [l3l and Rephaeli [l3| calcu- 
lated the photon frequency redistribution function in the 
electron rest frame usin g t he scattering probability de- 
rived by Chandrasekhar 3l|, which is called as Wright's 
method in the present paper. Another approach is the 
relativistic generalization of the Kompaneets equation^, 
where the relativistically covariant Boltzmann coUisional 
equation is solved for the photon distribution function. 
This approach was used by Challinor and Lasenby[l3| 
and Itoh, Kohyama, and Nozawa J^], which is called the 
covariant formalism in the present paper. Although the 
two are very different approaches, the obtained results 
for the SZ effect agreed extremely well. This has been a 
longstanding puzzle in the field of the relativistic study 
of the SZ effect for the last ten years. 

Very recently, Shimon and Rephaeli [s^jBoehm and 
J. Lavalle[3^ and Nozawa and Kohyama[33| (denoted 
NK hereafter) discussed the equivalence between differ- 
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ent approaches. In particular, NK showed that Wright's 
method and the covariant formalism were indeed mathe- 
matically equivalent in the approximation of the Thom- 
son limit, which is fully valid for the CMB photon en- 
ergies. This explained the reason why the two different 
approaches produced same results for the SZ effect even 
in the relativistic energies for electrons. Thus, NK clari- 
fied the situation of the longstanding puzzle. 

The third method for the study of the SZ effect is the 
direct numerical integration of the rate equation of the 
photon spectral distortion function. The first-order cal- 
culation in powers of the optical depth r was done by 
Itoh, Kohyama, and NozawajlSj for r ^ 1. The exact 
calculation which includes the full-order terms in pow- 
ers of T was done by Dolgov et al.fsS*! for r ^ 1. The 
rate equation in the present formalism is expressed in 
the form that is suitable to the direct numerical calcula- 
tion. In the present paper, we solve the rate equations 
numerically, and show the exact numerical solutions. 

In the present paper, we explore the formal solu- 
tions for the rate equations in the following represen- 
tation forms: the multiple scattering representation 34j . 
operator representation|36j, and Fourier transform 
r;presentationi3. In p^Tticular, we show that these 
representation forms are identical. With the operator 
representation form, we derive a formal solution that is 
equivalent to the formal solution obtained by the Fokker- 
Planck expansion approximation. 

The present paper is organized as follows: Starting 
from the rate equations derived in the NK paper, we de- 
rive in Sec. II the formal solutions of the rate equations 
for the SZ effect. In Sec. Ill, the present formalism is 
extended to the case which includes the peculiar veloc- 
ity of the CG. In Sec. IV, we solve the rate equations 
numerically, and obtain exact numerical solutions for the 
thermal SZ effect and kinematical SZ effect. The numeri- 
cal solutions for the nonthermal electron distributions are 
also presented in Sec. IV. Finally, concluding remarks are 
given in Sec. V. 



II. FORMAL SOLUTIONS FOR THE THERMAL 
SUNYAEV-ZELDOVICH EFFECT 

A. Rate equations 

In the NK paper, it was shown that the covariant 
formalism[l^ and Wright's method[l0[ were mathemat- 
ically equivalent in the approximation 7w/m <C 1 which 
is fully valid for the CMB photon energies where 
7 = 1 / -^Z 1 — , (3 and m are the velocity and rest mass 
of the electron, respectively. The two formalisms were 
connected each other by the Lorentz transformations for 
the zenith angles. Readers may be referred to the NK 
paperfsl] for the details. In the present paper, we use 
the expression in Wright's method for the photon fre- 
quency redistribution function, which was derived in the 
NK paper. 



In the present section, we consider the case that both 
of the CG and observer are fixed to the CMB frame. As 
a reference system, we choose the system that is fixed 
to the CMB. (Three frames are identical in the present 
case.) Throughout this paper, we use the natural unit 
h — c = 1 , unless otherwise stated explicitly. 

The rate equations for the photon distribution function 
n{x) and the spectral intensity function I{x) were derived 
in the NK paper, where x — Lo/ksTcMB is the photon 
energy in units of the thermal energy of the CMB. We 
recall the results here to make the present paper more 
self-contained. They are given as follows 34| : 



dn{x) 

dr 
dljx) 

dr 



dsPi{s) [n(e^x) - n{x)\ , 
dsPi{s) [liC^x) - I{x)\ , 

i 

dr — neCTdt , 



(1) 

(2) 
(3) 



where I{x) — Iox^n{x), Iq = (ksTcMB)^ /2tt^ , is the 
electron number density, gt is the Thomson scattering 
cross section. In Eqs. ([T]) and ([2]), -Pi(s) is the redistri- 
bution function for photon of a frequency shift s, which 
is defined by = x' /x, 



(4) 



Pl(s)= / dPP^^''pe{E)P{s,P), 

= ^T^sn / — ^t^/(mo,mo), (5) 



(6) 



where Pe{E) = m?pe{E) is the electron distribution 
function of a momentum p which is normalized by 
/(J dpp^pe [E) = 1 . The total probability is 



/C30 
dsPi{s) = 1. 
-OO 



(7) 



Variables appearing in Eqs. (|4]) - ([6]) are summarized as 
follows: 



/3™„ = (l-e-l^l)/(l + e~''') 
-1 



for s < 



[1 -e-"(l + /?)]//? for s > 



r [1-6-^(1 -/3)]//3 fors<0 
^2^*^ - 1 for s > 



(8) 
(9) 

(10) 
(11) 



It should be noted that the following useful relations 

P(5,/3)e-3- = P(-s,/3), Pi(s)e-^^ = Pi(-s) (12) 
were used in deriving Eq. 
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B. Formal solutions 

In the present section, we explore the formal solutions 
of the rate equations of Eqs. (P) and ^ in the following 
representation forms: </ie multiple scattering representa- 
tion, operator representation, and Fourier transform rep- 
resentation. In particular, we show that these represen- 
tation forms are identical. 

First, it is familiar to express the formal solutions in 
the multiple scattering representation. They are given, 
for example, in the NK paper 3J| as follows: 



n{x) 
I{x) 



dsP{s)no{e^x) , 
dsP(s)/o(e""x) 



(13) 
(14) 



where naix) and /o(a;) are the initial functions at r^O. 
In Eqs. and P{s) is expressed by 



i^(5) = e--5]^P,(.), 
j=0 ^' 
S{s), 



/OO nOO 
dsiPi(si)-- - / ds,_iPi(s,_i) 
CO J —OO 




(15) 
(16) 

(17) 



where Pj (s) is the redistribution function for the multiple 
scattering of the j-th order. 

Next, let us derive the formal solutions in the opera- 
tor representation. This method has an advantage that 
not only the solution can be expressed in a concise form 
but also it is useful to study relations between different 
methods. We write the following useful operator identity 
introduced by Bernstein and Dodelson|36[: 



-XD 



D 



d_ 

dx ' 



(18) 
(19) 



where A is a number and ,f{x) is an arbitrary function of 



Note that the 



factor goes to 1 since the deriva- 



tive has nothing to act on. It is important to understand 
the physical meaning of the operator D. Inserting A = s 
and f{x) = X into Eq. p8|) . one finds 



e X 



(20) 



where the definition of s (e" = x' /x) was used. Therefore, 
the operator D is nothing but the shift operator for a; in 
the present (Wright's) formalism. 
Introducing an operator 0{±D) by 



0{±D) 



dsPi{s)e 



±sD 



(21) 



and inserting into Eqs. H]) and ([2|), one obtains 
dn{x) 



dr 
dljx) 
~d7~ 



^0{D)n{x), 
= 0{-D)I(x)., 



(22) 
(23) 



where Eqs. (O and (|18p were used in the derivation. 
Thus, the following formal solutions in the operator rep- 
resentation are obtained: 



n{x) 



noix) , 



I{x)^e^<^^-^)h{x). 



(24) 
(25) 

It should be noted that Eqs. (|22p and ((23l) require the 
following condition for 0{D): 

x^O{D) =0(-D)x^ . (26) 

is guaranteed to be valid by 



In the present case, Eq. 
Eq. (Ull). 

We now show that Eqs. and (pS)) are identical 
to Eqs. (Uni) and respectively. Let us start with 

Eq. (gl) as follows: 



n{x) = exp <^T 

oc 



dsPi{s)e 



sD 



1 



noix) , 



dsFi(s)e*^ no{x). (27) 



Applying the identity of Eq. (|T8|) . the j-th order term 
(j ^ 1) can be rewritten as 



dsPi{s)e/"j no{x) 
dsiPi(si)-- - / ds,Pi{sj)no{e''+-+''x), 

-OO J —OO 

/OO pOO / i~ j- \ 

dsiPi{si) ■ ■■ dsPi s - '^o(e'*x) , 

-OO J-oo \ ^^-^ J 

dsPj{s)no{e''x) . (28) 



Inserting Eq. ([25)1 into Eq. ([?7|) . one finally obtains 
Eq. Uni). Thus, it has been shown that Eq. is iden- 
tical to Eq. The equivalence between Eqs. ([U]) and 
(1251) can be also shown in a similar manner. 

Let us now derive the third formal solutions for the rate 
equations. The Fourier transform representation was in- 
troduced by Taylor and Wright [33], where the redistri- 
bution function was derived in their Eq. (28). Consider 
the Fourier transform for P{s) of Eq. ([T5[) as follows: 



P{k) 



dsP{s)e 



— iks 



dsPj (s)e 



-iks 



(29) 
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Using Eq. (fT7|) , the j-th order term in Eq. ([29)l is written 
by 

/OO /'OO /'OO 

dsPj{s)e-''"' = dsiPiisi) ■ ■ ■ dsjPi{sj) 
-OO J —OO J —OO 



-ik{si-\ \-Sj) 



dsPi [s)e 



-iks 



(30) 



Inserting Eq. (f30|) into Eq. (|29| . one finally obtains the 
Fourier transform, 



P(fc) = 



,r[A(fc)-l] 



dsPi (s)e 



— iks 



The inverse Fourier transform for P(fc) is 



P(s) 



1 

2^ 



(31) 
(32) 



(33) 



Inserting Eq. ([55)1 into Eq. (fTH]) . one obtains 



n(a;) 



1 

2^ 



= gr[Pi{zD)-l] 



J —OO 

dke^^P^^'''>-^^S{k - iD)no{x) , 
no{x) , 



(34) 



where the operator identity of Eq. (|18p was used in 
the derivation. Thus, the formal solution for n{x) in 
the Fourier transform representation has been obtained. 
Similarly, the formal solution for I{x) is also obtained as 
follows: 



/(x) = e-[^i(-*^)-il/o(a;). 



(35) 



Here, we show that Eqs. (p4|) and ((35|) are identical to 
Eqs. and respectively. In deriving Eq. ([M)) . we 
found that 5{k ~ iD) was valid. By putting 

kD = iD, (36) 

and inserting into Eq. ([2T|) . one obtains 

/OO 
dsPi(s)e=F''=°'* - 1, 
-OO 

= Pl(±fcl3)-l, 

= Pi(±«i?)-l, (37) 

where Eq. p2p was used in the derivation. Inserting 
Eq. dniD into Eqs. ^ and (HS]), one finally obtains 
Eqs. ([34]) and ([35]) . Thus, the equivalence between the 
operator representation and Fourier transform represen- 
tation has been shown. Hence, the three representation 
forms are identical. 



C. Fokker-Planck expansion approximation 

Using the operator representation method, we ex- 
plore the relation of the formal solutions between the 
Fokker-Planck expansion approximation and the present 
method. For simplicity, we restrict ourselves to the SZ 
effect in the present paper, however, the extension to the 
kinematical SZ effect is straightforward. 

In Itoh, Kohyama, and NozawafTsj (denoted by IKN 
hereafter), the rate equation for the photon distribution 
function was expressed by their Eqs. (2.10)-(2.21) in the 
Fokker-Planck expansion approximation. Keeping the 
leading-order terms, the rate equation is reexpressed as 
follows: 



dn{x) 



ax 



2 5' r 3 4 

ox'^ ox-'' ox"^ 



9 5' . 10 r ' 



dx'^ 
n{x) 



8 <9« , 
ax^ 



(38) 



where Ii,I2t - ■ ■ ho are defined by Eqs. (2.12)-(2.21) of 
the IKN paper. The explicit forms are 

15.3 15 4,135 5 
1023^3 _ 2505^4 ^ 30375 



h = 40e 



el 



42 
^ 5 

11 3 

- 105^^ 
22 . 

^9 = 

7 

11 
210 



M06|2 



18594 ^. 

12059 ^. 

16568 

7516 

1 ( 

105 



7098 4 62391 



128 



'e I 



5 "^^^ 


10 


14253 4 


614727 


10 




40 


124389 


f 


+ 10 ' 




355703 


f 


+ 80 ' 





>5 

e 1 



e 7 



ho 



|5 

e ; 



(39) 
(40) 
(41) 
(42) 
(43) 
(44) 
(45) 
(46) 
(47) 
(48) 



where x — lo IksTcMB and O^. ~ kBTe/mc^ . It 
should be noted as follows: In Eqs. (2.12)-(2.21) of the 
IKN paper, terms not appearing in Eqs. (|39|) - (|48)) are 
0{TcMB/Te) <C 1, therefore, they are safely neglected. 

We now rewrite the rate equation of Eq. (|38p in terms 
of the operator D defined by Eq. (fT5|) . One obtains as 
follows: 



dn{x) 
dr 



OiKN{D)n{x), 

5 



(49) 
(50) 



5 



where 



(51) 



(52) 



31 
40' 



D 



1039 
120 



33 
10' 



11 



43 
10' 



269 

"so" 



(53) 



56 56 



2959 
140 



5993 
420 



1011 . 



140 



605 64^7 16 o , , 
84 35 105 ' ^ ' 



-D 



4371 

"896" 
3273 

80" 

14 



81707 



4480 



30339 



560 



4551 
560" 



3443. 



240 



199 



421 
105' 



11 

210' 



(55) 



Thus, the formal solution in the Fokker-Planck expansion 
approximation is obtained by 



(56) 



where no{x) is the initial solution at r = 0. 

On the other hand, the formal solution in the operator 
representation was given by 



n{x) 



0{D) 



dsPi{s)e'' 



1. 



(57) 
(58) 



By expanding e^^ in Eq. ([55)) . one obtains the following 

oo 

- J2 '^nD" , 



(59) 



n=l 



where 



— oo 
oc 



dsPi{s)s 

— / dpfpeifi) f d^Io 
nl Jo J-1 

f (Aio,Aio) 



1 



1 



■274(1-/3mo)3^ 



log 



dfi'o 

.1-/3Mo 



1 - /3mo 



(60) 



In deriving Eq. the definition of Pi{s) and = 

(1 — (3fiQ)/{l — Piiq) were used. Our task is to calculate 



Eq. ([601) for n = 1 to 10, and compare with Eq. dSOJ). The 
calculation is performed by expanding all functions in 
Eq. ([60|) in powers of /? and keeping up to 0(/3^°) terms. 
In the calculation, we used the symbolic manipulation 
program Mathematica. Then, the integration of Eq. (j60[) 
was also done analytically. The result is expressed by the 
power-series of 9^ ■ Finally, one has 



OiD)=OiKNiD) + Oi0t). 



(61) 



Thus, we have obtained a formal solution that is equiv- 
alent to the solution in the Fokker-Planck expansion ap- 
proximation. Similarly, the formal solution for the spec- 
tral intensity function in the Fokker-Planck expansion 
approximation is 



W(x) =e^°--~(-^)/o(x), 



(62) 



where Io{x) is the initial solution at t=0. It should be 
noted that the operator Oikn{D) satisfies the condition 
of Eq. ([26]). 



III. FORMAL SOLUTIONS FOR THE 
KINEMATICAL SUNYAEV-ZELDOVICH 
EFFECT 

A. Rate equations 

Let us now consider the case that the CG is moving 
with a peculiar velocity (3c {=Vc/c) with respect to the 
CMB. As a reference system, we choose the system that 
is fixed to the CMB. The z axis is fixed to a line connect- 
ing the observer and the center of mass of the CG. (We 
assume that the observer is fixed to the CMB frame.) In 
the present paper we choose the positive direction of the 
z axis as the direction of the propagation of a photon 
from the observer to the cluster. 

The rate equations for the photon distribution function 
n{x) and the spectral intensity function I{x) were derived 
in the NK paper[34], where x = Lo/ksTcMB- We recall 
the results here to make the present paper more self- 
contained. They are given as follows: 



dn{x) 

dr 
dljx) 

dr 



dsPi{s, l3c,z) [n{e''x) - n{x)] , 

dsP,{s,pc..) [e-3^/(e^a;) - /(x)] 

dr — UeaTdt , 
Pi(s,/?c,.) -Pi(s)+/3c,.Pi,k(s), 



(63) 

(64) 

(65) 
(66) 



where (3c,z is the peculiar velocity of the CG parallel to 
the observer. It should be noted that 0(/3c,z) ^^nd higher- 
order contributions were neglected in deriving Eq. (|66p . 
because (3c, z <C 1 is satisfied for most of the CG. The typ- 
ical value is /3c,z=l/300 for Wc=l,000km/s. In Eq. 
Pi(s) is the redistribution function defined by Eq. 
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Similarly, in Eq. (|66|) Pi^k{s) is the redistribution func- 
tion for photon of a frequency shift s in the case of the 
non-zero peculiar velocity of the CG, 



One has 



(67) 



1 



(1 - /3mo)3 



/(Mo,Aio), (68) 



where / (/io, Mo) given by Eq. ([5]), and variables appear- 
ing in Eqs. ^ and ^ are given by Eqs. ^ - pT|) . 
It should be noted that Eq. ([68]) was derived for elec- 
trons in thermal equilibrium at a temperature Tg. The 
distribution function for electron of a velocity /3 is given 
by 



Pe{E) = 



1 



(69) 



where 0^ — ksTe/mc^ , and K2{z) is the modified Bessel 
function of the second kind. For the power-law distribu- 
tions, J /Be in Eq. should be replaced by a//3^ and 
a for the p-power distribution Pe{E) cx p^'^ and the E- 
power distribution Pe{E) cx respectively. 



B. Properties of Pi(s) and Pi,k{s) 

In the present subsection, we study the properties of 
the redistribution functions Pi{s) and Pi,k{s)- It is fa- 
miliar that Pi{s) satisfies Eq. ([7]). We show that Pi,k{s) 
also satisfies 



/oc 
dsPi^Kis) = 1. 
-OO 



(70) 



if electrons are in a thermal equilibrium state. Integrat- 
ing Eq. (|67p over s, one obtains the following: 
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/OO PCX) 
dsPi,K{s) = / dpfp,{E) 
-OO Jo 

1 f+\ /?Mo-/?^ 



27*7-1 ^"(l-/3/io)*7-i 



c^Mo/(a*o,Mo) ,(71) 



where p = p/m. In deriving Eq. (j71[) . the integral vari- 
ables f3 and s were replaced by p and /Iq, respectively. 
Equation (|7ip is further simplified by inserting Eq. 
and using the following relations: 



1.2 ^7 



/OO /-OO -1 

dsPl.Kis) = / dpfp,{E)-fi^ 
-OO Jq 'J 



where (3 — p/E, j = E, and E = E/m were used. In the 
second line of Eq. ([7^ . {f{p)}' denotes the derivative of 
f{p) by p. Applying the partial integration to Eq. d74|) . 
one finally obtains 



dsPi,K{s) 



1 



dpp e 



1 



(75) 



Thus, Eq. ([70|) was shown. 

It is worthwhile to mention as follows: As seen from 
the right-hand side of the first fine of Eq. Eq. (|74|) 
is proportional to the pressure of the electron distribu- 
tion, because the pressure of the electron is described in 
general by 



' 1 
dpp^Pe{E)-pv , 



(76) 



where v — (3c. For the thermal distribution, one has 
Pthermai = nf,mc6f,, which again leads to Eq. (j70p . For 
the cases of nonthermal distributions, in general, Eq. ((70|) 
is not satisfied. This will be discussed in Sec. IV. 

More importantly, the present finding suggests that 
one can interpret Pi^k{s) as the pressure distribution 
function for a frequency shift s. In the present paper, 
we call Pi,k{s) the pressure distribution function here- 
after. 

In Figs, la and lb, we plot the frequency redistribu- 
tion function Pi (s) and the pressure distribution function 
Pi,k{s)^ respectively, as a function of s. The fc^Tg de- 
pendences are illustrated. The solid curve, dash-dotted 
curve, and dotted curve correspond to ksTe = 5keV, 
20keV, and lOOkeV, respectively. In Fig. la, the familiar 
behavior of the redistribution function Pi{s) is shown, 
where the sharp peak at s=0 for /cBTe=5keV becomes 
flatter and wider for higher electron temperatures. On 
the other hand, Pi^k{s) is shown for the first time. It 
is seen from Fig. lb that Pi.k{s) is not positive definite, 
which is interpreted as the pressure distribution. It is 
seen that (s) also has a large ksTe dependence. Two 
sharp peaks of opposite signs are canceling each other for 
the case kBTe=5keV, whereas the peaks become fiatter 
and wider for higher electron temperatures. 



c?Aio/(Aio,Mo) = 1, 



1 f+\ /?Mo - 1.2 



(72) 
(73) 



C. Formal solutions 

Let us now derive the formal solutions of the rate equa- 
tions of Eqs. and ([51)) in the operator representation. 
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lb 




FIG. 1: Plotting of Pi(s) and Pi,k{s) as a function of s. 
Figs, la and lb are -Pi(s) and Pi,k(s), respectively. The 
solid curve, dash-dotted curve, and dotted curve correspond 
to fcsTe = 5keV, 20keV, and lOOkeV, respectively. 



First, introduce the operators 0{~D) and OK{D — n) by 



/oc 
dsPi{s)e-''^ - 1 
-OO 

poo 

OK{D-n)^ / dsPi,K{s) 



(77) 

, (78) 



where n is an integer. Applying Eq. ((77l) and Eq. ((78|) of 
71 = 3 to the spectral intensity function I{x), one obtains 



0{-D)I{x)^ dsPi(s) [/(e~"x) -/(x)] , (79) 

J —OO 

/OO 
dsPi,K(s) 
-OO 



(80) 



Note that, in the case of the thermal distribution, 
Eq. ([75]) is further simplified as 

/OO 
dsPi,K(s)e^(^-")-l, (81) 
-OO 

because Eq. ([7D|) is valid. The rate equation of Eq. (IM)) 
is rewritten as follows: 



a/(a;) 



[0(-i?)+/3e,.0/f(i?-3)]/(x). 



(82) 



Now we expand the solution I{x) in powers of (3c. z as 
follows: 



lix) = Isz{x) + (ic,zhsz{x) + 0{Pl,) , 



(83) 



where SZ and fcS'Z denote the SZ effect and kinematical 
SZ effect of the first-order in /3c, z, respectively. Inserting 
Eq. ((83|) into Eq. ([82|) . one obtains the rate equation for 
each-order in /3c, z. 



dlszjx) 

dr 
diksz jx] 

dr 



= 0{-D)Isz{x), 
■=0{-D)hsz{x)+g{x,T), 



where 



g{x,T)^OK{D~i)Isz{x), 

-OK(i?-3)e^°(-^)/o(a;). 



(84) 
(85) 

(86) 



Equation (|84p corresponds to the rate equation for the 
SZ effect, which is identical to Eq. ([23l) . Therefore, the 
formal solution is given by Eq. (j25p . On the other hand, 
Eq. (|85p corresponds to the rate equation for the kine- 
matical SZ effect. 

Let us now solve Eq. (jSSp with an initial condition 
IkSz{x)=Q at T=0. Introducing a new function u{x) by 



Iksz{x) 



u{x) , 



(87) 



and inserting it into Eq. (|85p . one has the equation for 
u{x) as follows: 



du{x) 



dr 



-tO(-D) 



g{x,T) 



The solution is 



u{x) = / dr'e-"'<^(-^)5(2;,T'), 
Jo 



(88) 



(89) 



where the initial condition ii{x)~Q at t=0 was used. 
Thus, one has 



Iksz{x)^ / dr'e(---')'^(-^)g(x,r'), 
Jo 

Jq 

xe''°^-^ho{x). (90) 
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Equation (190]) is further simplified since the operators 
0{—D) and Ok{D — 3) are commutable each other. Fi- 
nally, one obtains the following formal solution for the 
spectral intensity function of the kinematical SZ effect: 

hsz{x) = tOk{D - 3)e"°(-^)/o(a:) , 
^TOK{D-3)Iszix), 
^rg{x,T). (91) 

Similarly, the rate equation for the photon distribution 
function n{x) is expressed by 



dn{x) 
dr 



= [0{D)+l3caOK{D)]n{x). 



(92) 



The initial condition is n{x) — no{x) at t—0. Expanding 
n{x) in powers of f3c.z by 



n{x) = nsz{x) + Pc,znksz{x) + 0(/3L) , 



(93) 



one finally obtains the formal solutions for the SZ effect 
and kinematical SZ effect as follows: 

nszix) = e^°(^)no(x) , (94) 
Ukszix) =TOK{D)nszix). (95) 

In deriving Eqs. ([M)) and the initial conditions 

nsz{x) = no{x) and nksz{x)~0 at t~0 were used. Be- 
fore closing this section, it should be noted that Eqs. (PT|) 
and (|95p require the following relation for Ok- 

x^Ok{D) ^Ok{D -i)x^ , (96) 

which is satisfied by Eq. ([75]) . 

IV. NUMERICAL SOLUTIONS 

A. Derivation of numerical solutions 

In this section we show the derivation of the numerical 
solutions for the rate equations of Eqs. (|84p and ((85)) . 
We consider an ideal condition that the CG is infinitely 
large. The rate equations are expressed as follows: 

^^''^^''^ - r dsP.is) [Iszie-^x) - Iszix)] , (97) 
dsPi{s) [Iksz{e~'x) - Iksz{x)] 



dr 
dikszix) 



dr 

+ g{x,T), (98) 
where g{x^T) is defined by Eq. (|86p . The expHcit form is 

/oo 
dsPi,K{s) [e-^'Isz{e'x) - Isz{x)] . (99) 
-oo 



The initial conditions at t = are given by 

^3 



Iszix) = /o(a:) = Iq- 



1 ' 



(100) 
(101) 



Iksz{x) =0, 

g{x,0) = / dsPi,K{s) [e-^'Io{e'x) - Io{x)] , (102) 



where /„ = (kBTcMB? /2tt^ . 

Let us first solve Eq. ([97)) . Introducing Isz{x,t) by 



Iszix) = e '^iszix,T) , 



(103) 



and inserting it into Eq. (|97p. one obtains the following 
equation for Iszix, t): 



diszix,' 
dr 



dsPiis)iszie-'x,T). (104) 



Integrating Eq. ()104p between r and t + At (At <C 1), 
and applying the Runge-Kutta method of the fourth- 
order, one obtains the following difference equation: 

Iszix, T + At) = iszix, t) 

+a(Ar) / dsPiis)iszie~'x,T), il05) 



liAr) = Ar + ^(Ar)^ + l(Ar)3 + l{Ar)^ . (106) 



(i) Inserting r = into Eq. (jlOSp , one has the solution 

aX T — At, 



iszix. At) = Ioix)+aiAT) / dsPi (s)/o(e-^x) , (107) 



where the initial condition of Eq. ()100p was used. In 
Eq. ()107p . the right-hand side (RHS) can be calculated 
with the initial function of Eq. ()100p . (ii) Similarly, in- 
serting T — At into Eq. (jlOSp . one obtains the solution 
at T = 2At, 

iszix, 2At) = iszix. At) 

nOO 

+a(Ar) / dsPiis)iszie~'x,AT),il08) 



where the RHS can be calculated with Eq. p07p . (iii) Re- 
peating the same step n-times, one obtains the solution 
at r = uAt, 

iszix, hAt) = iszix, in - 1)At) 

+a(AT) / dsPi(s)/sz(e-"x,(n-l)Ar),(109) 



where the RHS can be calculated with the solution at r = 
in — 1)At. Thus, the numerical solution Iszix) at r = 
NAt has been obtained. We call this method the full- 
order numerical calculation in contrast to the first-order 
numerical calculation in the r-expansion assuming r <^ 
1, which was done by Itoh, Kohyama, and Nozawa 15]. 

Similarly, Eq. (|98)) is also solved numerically. Intro- 
ducing IkSzix,T) by 



hszix) = e '^ikszix,T) 



(110) 
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and inserting it into Eq. ((98)) . one obtains the following 
equation for Iksz{x,T): 



2a 



dIksz{x,T) 
dr 

9{x,t) 



dsPi{s)iksz{e~"x, t) + g{x, t) ,(111) 

) 

-isz{x,T)\ . (112) 



It should be noted that g{x, r) in Eq. (jllip is a known 
function, because Isz{x,t) in Eq. (|112p was already ob- 
tained by solving Eq. (|104p . Integrating Eq. piip be- 
tween T and T + At{At <^ 1), and applying the Runge- 
Kutta method of the fourth-order, one obtains the fol- 
lowing difference equation: 

hsz{x, T + At) = iksz[x, t) 

+a{AT)i( dsPi{s)iksz{e~'x,T) 

[^J — OO 

+ ~9{x,t)\ . (113) 



(i) Inserting t — into Eq. (|113|) . one obtains the 
solution at r = At, 



iksz{x,AT) ^a{AT)g{x,0), 



(114) 



where the initial conditions of Eqs. poip and (I102p were 
used, (ii) Similarly, inserting r — At into Eq. (I113P , one 
obtains the solution at r = 2Ar, 

ikszix,2AT) ^ ikszix,AT) 

+a{AT)l dsPi{s)iksz{e-'x,AT) 

K.J — OO 

+ g{x,AT)\ , (115) 



where the RHS is calculated with Eq. (|114p . (ui) Re- 
peating the same step n-times, one obtains the solution 
at T = uAt, 

iksz{x,nAT) ^ iksz{x,{n - 1)At) 

+a( At) I / dsPi {s)iksz (e^^x, (n - 1) At) 

+ ~g{x,{n~l)AT)\ ,(116) 



where the RHS can be calculated with the solution at 
T ~ in — 1)At. Thus, the full-order numerical solution 
Iksz{x) at T — uAt has been obtained. It should be 
noted that Iksz{x) and nksz{x) can be also calculated 
directly from Eqs. ([91]) and ((95|) . respectively. 



B. Results for thermal distribution 

In this section we present the results in the full-order 
numerical calculations of Isz{x) and Iksz{x) for the 
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FIG. 2: Plotting of spectral distortion functions for t = 0.01. 
Figs. 2a and 2b are AIsz{x) and Alkszix), respectively. The 
solid curve, dash-dotted curve, and dotted curve correspond 
to fcsTe = 5keV, 20keV, and lOOkeV, respectively. 



thermal electron distribution. First we define the spec- 
tral distortion functions as follows: 

AIszix)/{Tlo) = (Iszix) - loix)) /{tIo) , (117) 
AIkszix)/iTlo) = Ikszix)/{Tlo) , (118) 



where loix) is defined by Eq. (|100p . 

We plot the spectral distortion functions as a function 
of X in Fig. 2. In Figs. 2a and 2b, the fc^Te dependences 
for AIsz{x) and AIksz{x) are studied for the case of 
T = 0.01. The solid curve, dash-dotted curve, and dot- 
ted curve correspond to the full-order calculations for 
A:BTe=5keV, 20keV, and lOOkeV, respectively. For most 
of the CG, fcsTe <20keV is satisfied, however, high ksTe 
is still needed when analyzing high temperature clusters, 
for example, Hansen, Pastor, and Semikoz[3^ showed 
that fcsTe >100keV is needed in order to get 2a contraint 
on the temperature. As seen from Fig. 2a, the peaks of 
AIsz{x) become higher and broader as ksT^ grows. In 
Fig. 2a, we also plot the results of the first-order numeri- 
cal calculations by Itoh, Kohyama, and Nozawa[15]. The 
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FIG. 3: Plotting of spectral distortion functions for 
fcsre=20keV. Figs. 3a and 3b are AIsz{x) and AIksz{x), re- 
spectively. The solid curve, dash-dotted curve, dashed curve, 
and dotted curve correspond to r = 0.01, 1, 10 and 20, re- 
spectively. 



two curves are indistinguishable in the entire region of 
X. This result guarantees that the first-order numerical 
calculation is valid for r ^ 1, which is satisfied for most 
of the CG. Similarly, in Fig. 2b the peaks of AIksz{x) 
become lower and broader as fesTe grows. In Fig. 2b, we 
also plot the results of the first-order numerical calcula- 
tions, where the two curves are again indistinguishable 
in the entire region of x. This result guarantees that 
the first-order numerical calculation is valid for the kine- 
matical SZ effect of the CG which satisfies r ^ 1. It is 
needless to mention that the signal of the kinematical SZ 
effect should be multiplied by the factor (3c. z- 

In Figs. 3a and 3b, we show the r-dcpendence of 
the spectral distortion functions as a function of x for 
fcBT'e=20keV. The solid curve, dash-dotted curve, dashed 
curve, and dotted curve correspond to t=0.01, 1, 10, and 
20, respectively. For most of the CG, r <C 1 is safely 
satisfied. On the other hand, r 3> 1 cases are very im- 
portant for optically thick plasmas, for example, in the 



accretion disk of the neutron star X-ray binaries [39|, and 
references are therein. It can be seen from Figs. 3a and 
3b that the r-dependence of the redistribution functions 
is large, where the present full-order calculations are in- 
dispensable for the cases r > 1. It can be seen that the 
peaks of AIsz{x) and Al^szix) become lower and ex- 
tremely broader as r grows. The full-order calculation for 
Isz{x) was also done by Dolgov et al.^j. The present 
calculation agrees with their calculation up to a;=l,000. 

In the present calculation, we have checked the conver- 
gence of the numerical solution for the step-size At. It 
has been found that At = 10""^ or less is necessary to get 
a good convergence, where the error was less than 0.05% 
compared with other numerical calculation. It was shown 
by Dolgov et al[35| that the solution has to be solved up 
to large values of x for the cases t ^ 1. This can be also 
seen from Figs. 3a and 3b. In order to study the con- 
vergence, we have calculated AIsz{x) and AIksz{x) up 
to a;=50,000 for t=20. The obtained ratios to the peak 
value at a;=l,000 and x=10,000 for Alszix) are 5.75% 
and 0.15%, respectively. Similarly, the ratios to the peak 
value at a:=l,000 and a;=10,000 for Ahszix) are 10.3% 
and 0.36%, respectively. 

Here, one should mention the difference between the 
present formalism and that of Dolgov et al.[3^ in the 
limit T — > oo. The Boltzmann equation in Dolgov et 
al. has, as a stationary solution, the thermal equilibrium 
one of the temperature Tg. On the other hand, the rate 
equations of Eqs. H]) and ^ which are derived in the 
approximation {to' — uj)/Te ^ 1 do not have the equi- 
librium solutions. For the equilibrium state, for exam- 
ple, of fcsTe ~ 20keV, the above condition is no longer 
valid. Thus, these equations are not applicable to the 
limit T —^ oo. The equilibrium solutions are recovered in 
the present formalism by replacing Eqs. ([1]) and ^ with 
new solutions n{x,r]) and I{x,r]) in the limit t — )■ cxi as 
follows: 



dn[x, rf) 
dr 



dr 



-n{x,Tj)\, (119) 
-/(a;,7y)l, (120) 



where 77 = TcMs/Te- In deriving Eqs. pT^)) and 
n(x, rj) <ti 1 was assumed in the original Boltzmann equa- 
tion. In the equilibrium state, this condition is justified, 
because the chemical potential |^| » 1. It can be seen 
that these equations take into account of the detailed 
balance corrections to Eqs. ^ and ([2]). The stationary 
solutions are 



I{x,T]) = Iox^n{x,T]) , 



(121) 
(122) 

where fi = lnC(3) + Sin?/ is the chemical potential in units 
of ksTe- The cross over frequency Xc in this limit is given 
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by Xc = — At, where r] = 1.174 x 10 * and = —54.60 at 
A:BTe=20keV. 

Finally, we derive the formal solutions for the improved 
rate equations. Introducing n[x, rj) = e~^^h{x) and 
I{x,T]) = e"''^/(a;), Eqs. (fTl9| and (fT20)) are rewritten 
as follows: 

dh{x) 



4a 



dr 

dijx) 
dr 



oo 



dsPi{s) [n{e^x) — n{x)] , 
dsPi{s) \i{e-'x) - i{x) 



(123) 
(124) 



where h(x) and I{x) satisfy the same rate equations as 
Eqs. (HI) and respectively. Thus, the formal solutions 
are 



n{x,r]) 



-rix^TO{D) 



no{x) 



(125) 

/(x, 77) = e-""e"°(-^)/o(x) . (126) 
where no{x) and Io{x) are the initial solutions at r = 0. 

C. Results for nonthermal distributions 

In this section, we calculate the spectral distortion 
functions for the nonthermal electron distribution func- 
tions. The careful study of the SZ effect on the non- 
thermal distributions was done, for example, by Co- 
lafrancesco et al.jl^, [s^]- In the present paper, we study 
two types of single power-law distributions for illustra- 
tive purposes. 

(i) the p-power distribution: 



Pe{E) 



An 



(flp - 



Pi <P <P2 

elsewhere 



Pi 



'P2 



(127) 
(128) 



where Eq. (|127p is normahzed by dpp'^pe{p) = 1. It 
should be noted that the constant Op is related to the 
constant a of Colafrancesco et al.[2^, [30| by ap = a + 
2, where we choose the reported values |29| a=2.5 and 
^2=1000, and study the pi-dependence of the spectral 
distortion functions, 
(ii) the £^-power distribution: 

Pe{E) = 7""^ , < p < 00 . 

. r(^ 



A-zy - — 



- 3 



(129) 
(130) 



where > 3, and Eq. I|129p is also normalized by 
dpp^Pe{p) — 1- In the present paper, we fix the pa- 
rameter value a-y=4.5 for an illustrative purpose. 

Using Eqs. p27p and (|129p . one can calculate the nor- 
malization of the pressure distribution Pi x(s)- One ob- 
tains 

/•oo 



dsPi^Kis) = 



1, 



p— power distribution (i^i) 
E— power distribution ' ^ ' 
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FIG. 4: Plotting of Pi{s) and Pi,i<r(s) for the nonthermal 
distributions. Figs. 4a and 4b are Pi{s) and Pi,k{s), respec- 
tively. The solid curve, dash-dotted curve, and dashed curve 
correspond to the p-power distributions for pi = 0.1, 0.3, and 
1.0, respectively. The dotted curve is the i5-power distribu- 
tion. In Fig. 4b, the dotted curve is multiplied by a factor 10 
in order to be visible in the same figure. 



where Eq. (|70p was satisfied for the iJ-power distribution 
case. 

First, in Figs. 4a and 4b, we plot Pi{s) and Pi^k{s), 
respectively. The solid curve, dash-dotted curve, and 
dashed curve correspond to the p-power distributions for 
Pi = 0.1, 0.3, and 1.0, respectively. The dotted curve is 
the i?-power distribution. In Fig. 4b, the dotted curve is 
multiplied by a factor 10 in order to be visible in the same 
figure. It is seen from Fig. 4a that -Pi(s) for pi=0.1 has 
a sharp peak at s=0 similar to the case of kBTf,=5keV 
in Fig. la, and the peak becomes flatter and wider as 
pi grows. Moreover, the peak position is shifting toward 
the positive direction. As seen from Fig. 4b, the struc- 
ture of Pi^K^s) is similar to the thermal distribution case, 
however, the parameter dependence is quite large. 
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FIG. 5: Plotting of spectral distortion functions for the non- 
thermal distributions at r=0.01. Figs. 5a and 5b are AIsz{x) 
and /S.Iksz{x), respectively. The solid curve, dash-dotted 
curve, and dashed curve correspond to the p-power distri- 
butions for pi = 0.1, 0.3, and 1.0, respectively. The dotted 
curve is the J5-power distribution. 



Similarly, in Figs. 5a and 5b, we plot l^Isz{x) and 
^Iksz{x) for r=0.01, respectively. The solid curve, dash- 
dotted curve, and dashed curve correspond to the p- 
power distributions for pi — 0.1, 0.3, and 1.0, respec- 
tively. The dotted curve is the _B-power distribution. The 
curves have long tails even for a small r values. More im- 
portantly, it is seen from Fig. 5b that Alksz shows a com- 
plex structure (two peaks) in the case of the i?-power dis- 
tribution. It is found that the similar structure appears 
for the p-power distributions for -pi > 0.8. This suggests 
that data for the kinematical SZ effect might give new 
limits to the parameter values for the nonthermal dis- 



tributions provided that the separation of the main SZ 
effect has been carefully done, for example, by taking into 
account for the non- isothermal profiles of the CG[40|. 

V. CONCLUDING REMARKS 

Starting from the rate equations for the photon distri- 
bution function which was derived in the NK paper, we 
derived the formal solutions for the SZ effect in three dif- 
ferent representation forms: the multiple scattering rep- 
resentation, operator representation, and Fourier trans- 
form representation. In particular, we showed that these 
representation forms were identical. By expanding the 
formal solution in the operator representation in powers 
of both the derivative operator D and electron velocity 
/3, we obtained the same formal solution derived in the 
Fokker-Planck expansion approximation. 

We also extended the present formalism to the kine- 
matical SZ effect. We obtained the formal solutions in 
the operator representation. Analytical properties for the 
photon frequency redistribution functions were studied. 
We found that the redistribution function which corre- 
sponds to the kinematical SZ effect can be interpreted as 
the pressure distribution of electrons. 

We solved the rate equation numerically, and obtained 
the exact numerical solutions for the thermal SZ effect 
and kinematical SZ effect. The exact numerical solutions 
include the full-order terms in powers of r, where r is 
the optical depth. We compared the present solutions 
with other calculations such as first-order calculation in 
T. As far as the clusters of galaxies are concerned (, where 
T ^ 1 is satisfied), the existing calculations are accurate 
enough to analyze the observational data. On the other 
hand, the exact calculation which includes the full-order 
terms in r is necessary for t > 1. 

Finally, we calculated the spectral intensity functions 
for the nonthermal electron distributions. The parameter 
dependences were studied for the p-power distribution. 
It was suggested that observation of the kinematical SZ 
effect might provide new limits to the parameter values 
for the nonthermal distributions. 
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